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. We present an effective potential and the corresponding modified Gross-Pitaevskii equation that 

account for the energy dependence of the two-body scattering amplitude through an effective-range 
expansion. For the ground state energy of a trapped condensate, the theory leads to what we call 
a shape-dependent confinement correction that improves agreements with diffusion Monte Carlo 
calculations. The theory illustrates, for relatively strong confinement and/or high density, how the 
' shape dependence on atom-atom interaction can come into play in a many-atom quantum system. 
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(L) . I. INTRODUCTION 

^ ^ The Gross-Pitaevskii equation (GPE) has in recent years played a central role in our understanding of gaseous 

d [ Bose-Einstein condensates (BEG) 0, • Its astonishing success has to do with the fact that in most experiments to 

■ date, both the gaseous parameter, na^ (here n is the number density and a is the scattering length), and the parameter 
I ; [ a/ Lq, which characterizes the strength of confinement (with Lq being a measure of the size of the trap), have been 

. very small. Recently, a new breed of experiments have emerged that take advantage of Feshbach resonances in cold 
^ atomic collision Q to make the scattering length a tunable through a magnetic field. Such experiments promise to take 
I us into parameter regimes in which the validity of GPE and its underlying assumptions, such as the shape-independent 
' ^ pseudopotential approximation have to be carefully examined. 

, Most of existing theories going beyond the GPE USB have focused on quantum fluctuation that is responsible 
for the lowest order correction in na^. But it has become clear that eventually, a better description of atom-atom 
interaction, beyond the pseudopotential approximation that is the basis of the GP theory, will be required. For 
homogeneous systems, this has been demonstrated in a recent work by Gowell et al. 0, which shows that different 
potentials having the same scattering length a can lead to vastly different ground state energy for na^ of the order 

,^ ^ of 0.05 or greater. For inhomogeneous systems, recent works on two atoms in a trap have shown that the shape- 

' independent approximation becomes less valid under strong confinement llOl . 

, Neither of these results are surprising. The shape-independent approximation assumes the two-body scattering 

■ amplitude to be a constant over the momentum (or energy) range of interest. For larger na^ or stronger confinement, 
[ a greater range of momentum states become involved and this assumption becomes less applicable. This point has 

■ already been illustrated for two- atom systems under strong confinement, where it has been shown that inclusion of 
the energy-dependence of scattering amplitude leads to much improved results for the energy levels of the system 
[3, Hoi . The goal of this article is to show how a similar idea can be generalized to a N-body system. 

In light of the success of the GP theory HQ, we focus here on a perturbative approach around the GPE. An effective 
potential is proposed that incorporates the energy dependence of the two-body scattering amplitude through an 
• effective-range expansion. A corresponding modified GP equation is derived. By taking advantage of the relationship 
f5 ' between effective range and scattering length as implied by the quantum-defect theory (QDT) 12, 13], no new inputs 
O from scattering calculations are required other than the standard scattering length and the Cq coefficient. For the 
ground state energy of an inhomogeneous system, our numerical results show that the modified equation leads to better 
agreements with diffusion Monte Garlo (DMG) results than theories that consider only the quantum fluctuation 
correction. More importantly, the theory introduces the concept of shape-dependent confinement correction, an origin 
of shape-dependence on the interaction potential that is absent in a homogeneous system. 
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II. MODIFIED GP EQUATION 

In a second-quantization formulation, a system of N bosons in an external trapping potential 14x1 (r) is described 
by a Hamiltonian 

J rfr$t(r) _|!.v2 + K=xt(r) $(r) + ^ J dn J dr2 $t(ri) $t(i.2)T/(ri _ r2)$(r2) $(ri) . (1) 
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Here V is the interaction between particles, ^'(r) and ^''''(r) are the bosonic annihilation and creation operators that 
satisfy the commutation relations 

[$t(ri), $t(i.2)] ^ [$(1.1), $(r2)] = 0, (2) 
[^{n),d\r2)]^5{r,-V2). (3) 
Taking the inter-particle potential V to be that of a mean field 

^MF = gS{ri - r2), (4) 

where 

9 = — — ' (5) 



and ignoring quantum fluctuation, one arrives at the usual GP equation 0, Q 

Hr,t). (6) 
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Here $ = ('5) and is nomalized by / dr|<i>p — N , with N being the total number of particles. 

From a two-body scattering point of view, the mean- field potential, given by Eq. is such that it gives, in the first 
Born approximation, the correct foward scattering amplitude at zero energy. It is a shape-independent approximation 
that ignores completely the energy-dependence of the scattering amplitude. While this may be a good approximation 
for small values of na'^ and for weak confinement, it becomes less valid for strong confinement and/or greater na^ . This 
is true even at zero temperature, as confinement and/or quantum depletion invariably involve nonzero momentum 
states. 

For a better treatment of atomic interaction, while preserving much of the structure of the GP theory, an effective 
interaction is proposed here that incorporates the energy dependence of the scattering amplitude through an effective- 
range expansion. Specifically, we look for a Vcs that gives, in the first Born approximation, the correct real part of 
the forward scattering amplitude, 

' dv e-^'^-'yeff 6+'"-^ = Re[/(fc, 9 = 0)1 (7) 



to the order of fc^ (the first order in energy). This requirement basically ensures a better optical potential for an atom 
moving in the medium of others p3 | . 

For the range of energies of interest here, only the s wave scattering is important. In this case, the scattering 
amplitude can be written in terms of the s wave phase shift (5o(fc) as 

■^^ ' ^ cot(5o(fc) — i ^ ^ 

From the standard effective-range expansion 

kcotSo{k) = + l-Vek'^ + (9) 
a 2 

where re is the effective range, we have, to the order of k'^, 

Re[fik,0)] = -a + a^a-^re)k\ (10) 

It is straightforward to verify that the Vcs satisfying Eq. Q can be written as 

V^cff = VMF(ri - r2) + Vmod (ri - r2), (11) 

where 

Knod = f [S{ri - r2)V2^_r, + ^l-rA^l - Ts)] , (12) 
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'iTih^a?{a - hr, 



(13) 



This potential is consistent with the pseudopotential expansion of Huang and Yang fl^l , except that it is Hermitian 
and is represented in the coordinate space. 

With this effective interaction, the many-body Hamiltonian becomes 



H — HmF + Hn^od 



Zm 



$(r)+i J dri J dr2$t(ri)$t(i.2)]4jp (ri - r2)$(r2) $(ri) 
+ ljdrij $t(ri) $t(r2)i>,„d (n - r2)$(r2) $(ri). (14) 

The term associated with Vmod can be simphfied in the two-body centre-of-mass frame [r = ri — r2, R = (ri + r2)/2] , 
in which i?mod takes the form: 

iJmod- ^32 j dnj dr$t(R+^)$t(R„|)[j(r)v^+V^<5(r)] $(R-|) $(R+£). 
After integrating over r and some integrations by parts, i?mod can be written as 



(15) 



i?mod = \g2 I dR*^(R) (*t(R)^(R)) ^(R). 



(16) 



Note that the Laplacian operator operates only on the number-density operator 4'^(R)\I'(R). Simplification of Hmf 
follows the usual procedure, and leads finally to 



Zm 



j dR$'^(R 
+ dR$+(R) [5($t(R)$(R) 



^(R) 



■ig2V2 ($t(R)$(R)) 



(17) 



From this Hamiltonian, both the static and the dynamic properties of a many-atom Bose system can be studied. We 
focus here on the ground state energy of the system. 

Ignoring both quantum fluctuation and iJmod, Eq. (|17|l leads to the standard GP energy functional 



(18) 



which corresponds to a time-independent GP equation for the ground-state wave function with a chemical potential 



Zm 



$(r,i) =/i$(r,t). (19) 
Inclusion of the lowest-order term due to quantum fluctuation leads to an energy functional Il6| 

^^MGPI [*] = J 



Zm 
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(20) 
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where the extra term due to quantum fluctuation was first derived by Lee, Huang and Yan^ 
the LHY term. This functional corresponds to a modified GP equation that we call MGPI 
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V^ + V,^t{r)+gm'+g,m' 



$(r,t) =/i$(r,t), 



and will be called 



(21) 



where 



91=9 



32 



,3/2 



(22) 



Inclusion of both quantum fluctuation and i?mod to the lowest order leads to an energy functional 



MGPII 



dr 



(23) 



The corresponding modified GP equation, which we call MGPII, is 

^S/' + V,Ar)+9m'+9im' 



2m 



-fv^(l 



<i>l 



$(r,i) =A^$(r,i). 



(24) 



MGPII incorporates the shape dependence on the interaction potential through the effective range r^. This addi- 
tional parameter does not, however, add much extra complexity to the theory, as and a are generally related. For 
a hard sphere potential, we have the well-know relation 



For atoms with van der Waals interaction. 



re = (2/3)a. 
can be determined from a by [l 

3XeJ {a/f^Q 



1 



(«//36)]'}, 



(25) 



(26) 



where (3^ — {mCQ/h^Y^'^ '"^ length scale associated with the van der Waals interaction, and Xe = [r(l/4)]"^ / (27r). 
Similar relations for other long-range potentials should also exist, an assertion that can be deduced from a more 
general QDT consideration 13]. 

The physical meaning of the correction due to -ffmod can be better understood by estimating the ordcr-of-magnitude 
of each term in Eq. (|23l) or (|24|) . Consider a system of N bosons confined to a length scale Lq. The ordcr-of-magnitude 
for the number density is n ^ N/{Lq)^. Let cmf ^ na'^(ft^/2m)(l/a)^ to represent the energy per particle due to the 
mean field. It is easy to show that the LHY quantum fiuctuation term is of the order of eLHY ~ £MF("-a'^)^^^: while 
the i?mod term is the order of e,jiod ~ eMF(a/io)^- It is thus clear that the two corrections are of different origin. The 
Hnwd correction has an order-of-magnitude that is determined primarily by the strength of confinement. It depends 
also on the shape of the interaction potential through g2 and rg. We will call this correction the shape-dependent 
confinement correction. It is a new source of shape dependence that is absent in a homogeneous system (Lq — > oo with 
n being fixed). In this regard, it differs substantially from the shape dependence that can come at higher densities 
from two-body correlation at shorter length scales ja, 113 ■ 

A further comparison of the magnitudes of the LHY and the shape-dependent confinement corrections reveals 
other interesting features. From the estimates above, we have emod/eLHY ~ {a/ LqY / [na^Y^"^ ~ {a/ Lq)^/"^ /N^/'^ . 
This means, for instance, that relative to the LHY correction, the shape-dependent confinement correction is more 
important for small N than for large TV, a point that will be further demonstrated numerically. 

We have chosen, up to this point, a formulation that most closely parallels the standard GP formulation 
a small particle number iV, a number-conserving Hartree-Fock formulation is in fact more appropriate llE 
corresponding results are easily derived. For example, at the level of MGPII, we have 

|V0|Vyext(r)|0|' 




MGPII 



[(f)]^ N j dr 
+ l9iN- 
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l9.iN- 
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:92iN^ 



1)|</.|'V2 



(l^r) 



(27) 
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FIG. 1: Relative difFerences between ground state energies predicted by the GP theory and by diffusion Monte Carlo calculations 
0- The results are for A'' bosons in a spherically symmetric trap of size Lq — \J h/2mu. Here a is the scattering length. 
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+ Kxt(r) + 9{N - 1) \^\^ + g^{N - if'^ |0|' 
(iV-l)V2(|0|2)]0(r,t)=/.0(r,i), 



(28) 



where is normalized according to Jdrj^p = 1. It is the Hartree-Fock formulation that we actually implement in 
our numerical calculations, but the difference is small for large N (between N and — 1). Setting (72, or both §2 and 
f/i, to zero in Eqs. H27|l and (|28|l leads to Hartree-Fock results in MGPI and GP, respectively. 



III. NUMERICAL RESULTS 



We present here our numerical results for the ground state energy of a BEG in different approximations. They 
are obtained using an imaginary time evolution method with ADI algorithm T^. The code has been tested by 
reproducing the GPE results of Ref . 2Q\ . 



A. Results for hard spheres 

A many-body system of hard spheres is an ideal test-ground for various theoretical models, as accurate results for 
its ground-state energy are available from diffusion Monte Garlo (DMG) calculations for both homogeneous and 
inhomogeneous cases 

Figure 1 shows the relative differences of the energy functional between GP and DMG for N bosons in a spherical 
harmonic trap characterized by a length scale Lo = \/ h/2mLo. It is clear from these results that the ground state 
energy of the system can differ substantially from those predicated by the GP theory, for either strong confinement 
(large a/Lo) and/or high density (large N). 

As shown in Ref. 7], most of these differences, at least for the range of parameters considered, can be accounted for 
by the LHY quantum fluctuation correction at a level corresponding to MGPI. This is illustrated in Figure 2(a). 
In the same figure, it is also shown that a consistently better agreement with DMG is achieved by MGPII, which 
includes the shape-dependent confinement correction. This means that of the remaining difference between MGPI 
and DMG, at least part of it is shape-dependent and requires a better description of atom-atom interaction. Note 
that in the case of hard spheres, (72 is always positive. The sign of the shape-dependent confinement correction is then 
determined by the expectation value of (l*^!^) ' "^hich is negative for the ground state. 

More results of MGPII are presented in Figure 2(b). It shows that at the level of MGPII, the relative difference 
from DMG has been reduced from up to 20% for the GP theory (see Figure 1) to less than 2% for the range of 
parameters considered. 
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FIG. 2: Relative differences between ground state energies predicted by modified GP tlieories and by diffusion Monte Carlo 
calculations 0- (a) A sample comparison of MGPI and MGPII, both of which lead to much better results than the GP theory, 
(b) More results of MGPII. Note that the t/-scale of these figures are about a factor of ten smaller than that of Figure 1. 
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FIG. 3: The universal relationship between r^/ and a/ fi& for atoms with long-range van der Waals interaction [l^ . 



B. Results for atoms with van der Waals interaction 

With the introduction of -ffmod, the ground state energy becomes dependent upon the shape of the long-range 
potential through the relationship between Ve and a. For atoms with van der Waals interaction, r,, and a are related 
by Eq. H26() . Figure 3 is a graphic illustration of this relation. Note that diverges for a ~ 0, corresponding to the 
fact that in this particular case, the energy dependence of the scattering amplitude around zero energy cannot be 
described by an effective-range expansion |22l |. Also note that unlike the case of hard spheres, the g2 for a van der 
Waals potential may become negative for sufficiently small a/ (3q. Consequently, the shape-dependent confinement 
correction to the ground state energy may in principle be positive. [Both points suggest that something interesting 
happens around a = that may deserve a separate investigation.] 

Figure 4(a) gives a comparison of ground state energies predicted by MGPII and by the GP theory for a realistic 
experimental configuration. Specifically, it is for a *^Rb condensate in a cylindrical trap with an aspect ratio of 
£ = 6.8/17.5 23J. A Cq value of 4698 a.u. is used (2^. It shows that deviations from the GP theory can become 
substantial for either high density and/or strong confinement. With the tunability of the scattering length via a 
Feshbach resonance |2^, such many-body effects beyond the GP theory may soon become observable. Figure 4(b) 
shows more explicitly the shape-dependent confinement correction. Its magnitude increases with a/Loj as expected. 
For a fixed a/Lo, its relative contribution is more significant for small N than for large N , consistent with our earlier 
discussion. The corrections are mostly negative (except for very small a) since 52 is greater than zero for most of 
the data points shown. [Because Lg = 34850 a.u. is much greater than — 164.2 a.u., most of the data points 
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FIG. 4: (a) Relative differences between ground state energies predicted by MGPII and by tlie GP tlieory for a *^Rb condensate 
in a cylindrical trap of aspect ratio e = 6.8/17.5. Here Lo is defined by the transverse frequency: Lo = y^y2ma7L. (b) Relative 
contribution from the shape-dependent confinement correction. 

correspond to a/P^ ^ 1. The is then approximately a constant « (37r)^^[r(l/4)]^/36, and the corresponding g2^s 
are positive.] 

IV. CONCLUSION 

In conclusion, an effective interaction and the corresponding modified GP equation have been proposed that take 
into account the energy dependence of the two-body scattering amplitude through an effective-range expansion. The 
resulting theory, called MGPII, leads to better agreements with diffusion Monte Carlo calculations [3 than either the 
GP theory, or MGPI that considers only the quantum fluctuation correction. The theory expands, considerably, the 
parameter space (specified by na^ and a/Lo), in which a GP type of formulation can be applied. It introduces the 
concept of shape-dependent confinement correction and shows that for a fixed confinement, its relative contribution 
is more significant for small N than for large N. 

Finally, we point out that both MGPI and MGPII are perturbative expansions around the GP theory. They will 
eventually fail for sufficiently large na^ and/or sufficiently strong confinement. For high density, shape dependence 
of a different origin will eventually emerge from two-body correlations at short length scales [3, [131 . For strong 
confinement, the local density approximation, implicit in both MGPI and MGPII, will eventually fail. Exploration 
into those regimes will require theories that differ much more substantially from the GP formulation. 
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